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Abstract 

Background: Cadmium (Cd) exposure and sulfate limitation induce root sulfate uptake to meet the metabolic 
demand for reduced sulfur. Although these responses are well studied, some aspects are still an object of debate, 
since little is known about the molecular mechanisms by which changes in sulfate availability and sulfur metabolic 
demand are perceived and transduced into changes in the expression of the high-affinity sulfate transporters of the 
roots. The analysis of the natural variation occurring in species with complex and highly redundant genome could 
provide precious information to better understand the topic, because of the possible retention of mutations in the 
sulfate transporter genes. 

Results: The analysis of plant sulfur nutritional status and root sulfate uptake performed on plants of Brassica juncea - a 
naturally occurring allotetraploid species - grown either under Cd exposure or sulfate limitation showed that both 
these conditions increased root sulfate uptake capacity but they caused quite dissimilar nutritional states, as indicated 
by changes in the levels of nonprotein thiols, glutathione and sulfate of both roots and shoots. Such behaviors were 
related to the general accumulation of the transcripts of the transporters involved in root sulfate uptake {BjSultrl;! and 
BjSultrl;2). However, a deeper analysis of the expression patterns of three redundant, fully functional, and simultaneously 
expressed Sultrl;2 forms {BjSultrl ;2a, Bj5ultrl;2b, BjSultrl;2c) revealed that sulfate limitation induced the expression of all 
the variants, whilst BjSultrl ;2b and BjSultrl;2c only seemed to have the capacity to respond to Cd. 

Conclusions: A novel method to estimate the apparent k^ for sulfate, avoiding the use of radiotracers, revealed that 
BjSultrl;! and BjSultrl ;2a/b/c are fully functional high-affinity sulfate transporters. The different behavior of the three 
BjSultrl ;2 variants following Cd exposure or sulfate limitation suggests the existence of at least two distinct signal 
transduction pathways controlling root sulfate uptake in dissimilar nutritional and metabolic states. 
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Background 

Sulfur is an essential element for all living organisms, since 
it is found in a broad variety of biological compounds play- 
ing pivotal roles in a number of metabolic processes [1]. In 
contrast to animals, which have a dietary requirement for 
some organic sulfur compounds, plants have metabolic 



* Correspondence: fabio. nocitoOunimi. it 

'Dipartimento di Scienze Agrarie e Ambientali - Produzione, Territorio, 
Agroenergia, Universita degli Studi di Milano, 20133 Milano, Italy 
Full list of author information is available at the end of the article 

(3 Bion/led Central 



pathways that allow them to assimilate inorganic sulfur 
into organic sulfur compounds through a cascade of well 
characterized enzymatic steps. For this reason plant sul- 
fur assimilatory pathways are considered to be the main 
sources of organic sulfur compounds for animal and hu- 
man diets [2]. 

The main sulfur source for plants is the sulfate ion of 
the soil solution available in the rhizosphere [3,4], which is 
taken up through specific root plasma membrane high- 
affinity sulfate transporters. Once inside the plant, sulfate 
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is allocated to different sinks, and undergoes intracellular 
channeling to chloroplast and vacuole, where it is assimi- 
lated into organic sulfur compounds or compartmental- 
ized as sulfur store, respectively [2]. The main pathway 
of sulfate assimilation in plants involves the adenylation 
of the anion and its stepwise reduction to sulfite and then 
sulfide which is finally incorporated via 0-acetylserine 
(OAS) into cysteine (Cys), a key intermediate from which 
the essential amino acid methionine (Met), the tripeptide 
glutathione (GSH), and most sulfur containing compounds 
are synthesized [2,5]. 

Considering the central role of Cys in sulfur metabolism, 
it appears evident that both sulfate uptake and the reduc- 
tive assimilation pathway have to be finely modulated to 
meet the metabolic demand for sulfur arising from Cys 
consuming activities, which largely contribute to define 
the total sulfur requirement of plants. Such a demand may 
consistently vary under the different environmental condi- 
tions that plants may experience during their growth. For 
instance, biotic and abiotic stresses may increase the meta- 
bolic demand for some Cys derived compounds, causing 
an increase in the activity of the sulfate assimilatory path- 
way [6]. An example of this has been largely described in 
plants exposed to cadmium (Cd) in which the activation of 
a wide range of adaptive responses involving GSH con- 
suming activities may increase the demand for sulfate, sul- 
fur metabolites and carbon skeletons [7-10]. Indeed, GSH 
not only acts as an antioxidant in mitigating Cd-induced 
oxidative stress, but also represents the key intermediate 
for the synthesis of phytochelatins (PCs), a class of Cys-rich 
heavy metal-binding peptides involved in buffering cytosolic 
metal-ion concentration [11]. The large amount of PCs 
produced by Cd stressed plants represents an additional 
sink for reduced sulfur which, by increasing the metabolic 
request for both Cys and GSH, generates a typical demand- 
driven coordinated transcriptional regulation of genes 
involved in sulfate uptake, sulfate assimilation and GSH 
biosynthesis. Such a response is thought to be essential 
to satisfy two contrasting needs arising from Cd stress: 
i) maintaining cell GSH homeostasis; ii) detoxifying heavy 
metals by means of GSH-consuming activities. A similar ac- 
tivation has been described under sulfate limitation [12-14], 
although in this condition plant sulfur needs to sustain the 
growth do not vary: the induction of sulfate transporters 
and enzymes along the assimilatory pathway reflects some 
difficulties in maintaining both an adequate rate of Cys bio- 
synthesis and sulfur-containing compound homeostasis. 

Sulfate transport activations under Cd stress and sul- 
fate limitation have been shown to be mainly controlled 
at transcriptional level and have been often indicated as 
resulting from the same, although controversial, nutri- 
tional signals [8,9,15]. In the current model of transcrip- 
tional regulation, some intermediates along the pathway of 
sulfate assimilation and GSH biosynthesis act as negative 



or positive signals in modulating the expression of sulfate 
transporters. Adequate levels of reduced sulfur compounds, 
such as Cys and GSH, would repress gene expression 
through a negative feedback loop preventing excessive sul- 
fate uptake and reduction; vice versa a contraction of GSH 
pools would de-repress gene transcription allowing sulfate 
to enter the pathway. A second regulatory loop, involving 
OAS as a key intermediate, should act in promoting gene 
de-repression when nitrogen and carbon supply exceeds 
sulfur availabOity within the cells. In this condition, since 
sulfide avaHabUity is not enough for Cys biosynthesis, OAS 
accumulates and partially overrides the negative feedback 
provided by GSH on gene transcription [16]. Such a revers- 
ible regulation allows the system to adjust sulfate uptake 
to the nutritional status of the plant, and agrees with the 
concept of demand-driven regulation of sulfate uptake 
and metabolism [12]. 

Comparative studies clearly show that both sulfate 
deprivation and Cd stress produce a contraction in the 
GSH pools and a positive change in the OAS levels, which 
in turn may induce the accumulation of high-affinity sul- 
fate transporter mRNAs, allowing sulfate to enter the cells 
[15]. However, some aspects of this picture need to be fur- 
ther investigated, since the relationships existing between 
the accumulation of sulfate transporter mRNAs and the 
levels of the signal-intermediates do not always appear to 
be evident [9,17]. Moreover, Rouached and co-workers 
[15] clearly showed that the expression of the Arabidopsis 
Sultrl;! and Sultrl;2 - two high-affinity sulfate transporter 
genes - is not regulated in complete agreement with the 
current model, and they proposed the existence of distinct 
signaling pathways controlling sulfate uptake under differ- 
ent sulfur nutritional status. Finally, whether cellular con- 
tents of sulfate, sulfide, OAS, Cys and GSH are the true 
primary signals for controlling sulfate uptake and reduc- 
tion or rather act indirectly is still a matter of investigation 
[14,18], since very little is known about the molecular 
mechanisms involved in the nutritional signal perception 
and transduction [2,19,20]. Thus the need for additional 
efforts and integrated experimental approaches appears 
particularly evident to unveil this picture. The analysis of 
the natural variation occurring in species with redundant 
genomes could provide precious information about the 
molecular mechanisms controlling sulfate uptake, since 
the presence of redundant genes may have led to the accu- 
mulation of mutations which otherwise would have been 
eliminated by natural selection. From this point of view 
the species belonging to the Brassica genus could be very 
useful, since several lines of evidence suggest that the 
genomes of the three diploid Brassica species {B. rapa, 
B. oleracea and B. nigra) are composed of three rear- 
ranged variants of an ancestral genome - structurally simi- 
lar to that of Arabidopsis thaliana - and descended from 
a common mesohexaploid ancestor [21-23]. Moreover the 
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level of complexity may be further increased by considering 
the allopolyploid Bmssica species in which two distinct 
Bmssica genomes cohabit [24], increasing the probability of 
evolving novel gene interactions through the processes of 
sub-functionalization and/or neo-functionalization of para- 
logs [25,26]. 

In this work we present and discuss some evidence 
toward the existence of multiple transduction pathways 
controlling sulfate uptake under Cd stress and sulfate 
limitation in Bmssica juncea (AABB, n = 18), a natural oc- 
curring allotetraploid species formed through hybridization 
between B. rapa (AA, n = 10) and B. nigra (BB, n = 8), as 
described by the "triangle of U" [24]. 

Methods 

Plant material, growth conditions, and 
experimental design 

Brassica juncea L. Czern & Coss (Lodi selection) seeds 
were sown on filter paper saturated with distilled water 
and incubated at 26°C in the dark. Three days after sowing, 
seedlings selected for uniform growth were transplanted 
into 5 L plastic tanks (6 seedlings per tank) containing an 
aerated complete nutrient solution [500 i^M NH4H2PO4, 
3 mM KNO3, 2 mM Ca(N03)2, 1 mM MgS04, 25 
Fe-tartrate, 46 |iM H3BO3, 9 ^iM MnClz, 0.8 ZnCla, 
0.3 nM CuCla, 0.1 nM (NH4)6Mo7024, pH 6.5] and kept 
for 14 days (pre-growing period) in a growth chamber 
maintained at 26°C and 80% relative humidity, with a 
16-h light period. For Cd treatments, plants were grown 
for an additional 8 days (acclimation period) in a 5-fold 
diluted (not for micronutrients) nutrient solution (accli- 
mation solution) and then exposed to different Cd con- 
centrations (0, 10, and 25 CdQa) for 48 h. For sulfate 
limitation treatments, at the end of the pre-growing period 
plants were grown for 10 days in the acclimation solu- 
tion containing different sulfate concentrations (200, 
50 or 10 i^M); in the cases of the lowest sulfate con- 
centrations, MgCl2 was added to maintain the same 
concentration of magnesium. In both cases the growth 
chamber parameters were the same as described be- 
fore, and all hydroponic solutions were renewed twice 
a week to minimize nutrient depletion. At the end of the 
experimental periods, plants were immediately used for the 
in vivo experiments or harvested to be further analyzed. In 
this case roots were washed for 10 min in ice-cold 5 mM 
CaCl2 solution to displace extracellular Cd [27], rinsed in 
distilled water and gently blotted with paper towels; shoots 
were separated from roots and the tissues were fi-ozen in 
liquid N2 and stored at -80°C. 

RNA extraction and cDNA cloning 

BjSultrl;! and BjSultrl;2 partial cDNAs were amplified 
by RT-PCR from Brassica juncea mRNA isolated from 
roots. Total RNA was extracted from roots of sulfur- 



starved plants using TRIzol reagent (LifeTechnologies), 
poly A* mRNA was isolated using the Oligotex mRNA 
Spin-Column system (QIAGEN), and first-strand cDNA 
synthesis was carried out using the SuperScriptlll first- 
strand synthesis system for RT-PCR (LifeTechnologies) 
according to the manufacturer's instructions. Degenerate 
primers BjSultrljldegdir (5'-ACGGAGGAGGGTCCGRTG 
CAA-3'), BjSultrl;ldegrev (5'-TTYGGGTCGATCACGGCC 
TGGCA-3'), BjSultrl;2degdir (5'-GTYTTCGATTGGGGRC 
GTAR-3'), and BjSultrl;2degrev (5'-RAGGAAGAGCAATG 
TCAAGAGA-3'), were designed based on highly con- 
served regions identified in sequences of sulfate trans- 
porter cDNAs of Brassica napus and Arabidopsis thaliana 
[for BjSultrl;!: BnSultrl;! (GenBank accession no. AJ41 
6460) and AtSultrl;! (TAIR accession no. At4g08620); for 
BjSultrl;2: BnSultrl;2 (GenBank accession no. AJ311388), 
and AtSultrl;2 (TAIR accession no. Atlg78000)]. 5'- and 
3 '-regions of the sulfate transporter cDNAs were iso- 
lated by 5 - and 3 -RACE approach using GeneRacer 
Kit (LifeTechnologies) according to the manufacturer's 
instructions. Finally the full coding regions were con- 
firmed by RT-PCR using sequence specific primers 
obtained from the 5 - and 3 -RACE fragments, and 
proofreading PJu-DNA polymerase (Promega). All PGR 
products were verified by sequencing after cloning into 
the pCR-Bluntll vector (LifeTechnologies), and sequence 
data were submitted to GenBank (accession no. JX896426, 
BjSultrl;!; JX896427, BjSultr!;2a; JX896428, BjSultr!;2b; 
JX896429, BjSultr!;2c). 

Sequence analyses were performed using ClustalW 
and neighbor-joining trees were generated using MEGA 
5.05 [28]. 

Gene expression analysis 

Semi-quantitative RT-PCR analyses of BjSultr!;! and 
BjSultr!;2 pool were performed on first-strand cDNA 
deriving from total RNA extracted from roots. PGR was 
carried out for 24 cycles, where cDNAs were exponen- 
tially amplified by 7^-DNA polymerase (Promega), using 
the following couples of primers: BjSultrl;ldir 5'-ACGG 
AGGAGGGTCCGATGCAA-3' and BjSultrl;l,ev 5'-TTC 
GGGTCGATCACGGCCTGGCA-3' (producing a 453 bp 
fragment), BjSultrl;2dir 5'-GGTTTTCGATTGGGGACG 
TA-3' and BiSultrl;2rev 5'-TGTCAAGAGAACAACGATT 
GAC-3' (producing 1046 bp overlapping fragments). cDNA 
loading was normalized using the BjTub 846 bp amplicon 
(accession no. JX896430), as an internal control, obtained 
with primers designed on conserved regions of beta tubulin 
Tub9 sequences of Arabidopsis thaliana (TAIR accession 
no. At4g20890) and Brassica napus (GenBank accession 
no. AF258790) as follow: Tubjir 5'-TGTTGTGAGGAAG 
GAAGCTGAG-3' and Tub,ev 5'-TCCTGTGTACCAATG 
AAGG-3 '. PGR products were separated in agarose gels 
and stained with SYBR Green I (LifeTechnologies); signals 
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were detected using a laser scanner (Typhoon 9200, GE 
Healthcare) with a 532 nm laser and a 526 nm filter. 

For semi-quantitative RT-PCR analyses of the three dif- 
ferent variants of BjSultrl;2, the entire ORFs were amplified 
with BiSultrl;2ATG 5'-ATGTCTGGGAGAGCTCATCCT 
G-3' and BjSultrl;2sTOP 5'-TCAGACCTCGTCGGAGAG 
TTTTG-3' primers (producing a fragment of 1968 bp for 
BjSultrl;2a and fragments of 1959 bp for BjSultrl;2b and 
BjSultrl;2c). PGR products were then digested with Clal 
endonuclease at 37°C for 3 h, and restriction products 
were separated in agarose gels. Signals were detected 
after staining as above described, and densitometrically 
analyzed using ImageJ 1.46 software [29]. 

All the expression analyses were performed using three 
independent cDNAs deriving from three independent ex- 
periments in which six plants were pooled for RNA extrac- 
tion. Each cDNA was amplified, digested, run on gel, and 
quantified three times (m = 9). 

Heterologous expression of sulfate transporters and 
kinetic analysis in yeast 

EcoRTended fragments, resulting from the amplification 
of BjSultrl;!, BjSultrl;2a/b/c, ZmSTl;l, and AtSultr2;l 
ORFs using appropriate primers (BjSultrl; Ikatg 5'-CA 
CTAGAATTCTAAAAAATGGCCAAGACTAATCCGC 
CGGA-3' and BjSultrl;lKSTop 5'-TGACCGAATTCTT 
ATGCTTGTTGCTCAGCCAAT-3', BjSultrl;2icATG 5'- 
CACTAGAATTCTAAAAAATGTCTGGGAGAGCTCAT 
CCTG-3' and BjSultrl;2KSTOP 5'- TGACCGAATTCTCA 
GACCTCGTCGGAGAGTTTTG-3', ZmSTl;lKATG 5 -C 
AGCGAATTCTAAAAAATGCCGCCGCGAACGGTGTC 
C-3' and ZmSTl;lKSTOP 5'-GCGCGAATTCTCAGACAT 
TATCGACCATCTTAGGAGC-3', and AtSultr2;licATG 5'- 
CAGCGAATTCTAAAAAATGAAAGAGAGAGATTCAG 
AGA-3' and AtSultr2;lKSTOP 5'-TGACCGAATTCTTAA 
ACTTTTAATCCAAAGCAAGCATCAA-3') including a 
consensus sequence for translation initiation in yeast 
[30], were subcloned in the EcoRI site of the yeast {Sac- 
charomyces cerevisiae) expression vector pESC-TRP 
(Stratagene) under the control of GALIO promoter. 
Chimeric and empty vectors were used to transform the 
yeast double sulfate transporter mutant CP154-7A (MATa 
his3 leu2 um3 ade2 trpl sull:LEU2 sul2:URA3) [31] using 
the standard lithium acetate method [32], and Trp* recom- 
binant yeast cells were selected. Complementation tests 
were performed as previously described [9] . 

For the growth analysis, recombinant yeast cells were 
grown - at 28°C in a synthetic Trp-free liquid medium 
containing yeast nitrogen base and required amino acids - 
up to reach a mid-log phase. Yeast cells were then washed 
twice with sterile distilled water and resuspended to a final 
absorbance of 0.1 ^eoo unit in the B minimal medium [33], 
supplemented with 40 |ig mL"^ adenine and 200 i^g mL"^ 
histidine to meet the auxotrophies of the strain, and 



containing different amounts of Na2S04 or 100 [iM DL- 
homocysteine (HCys) as sole sulfur sources. Yeasts were 
incubated at 28°C and growth was monitored by measur- 
ing the absorbance at 600 nm. At the end of the growing 
period, 30 mL of the yeast culture was harvested, washed 
twice in sterile distilled water, resuspended in 4 mL of boil- 
ing buffered ethanol (75% ethanol in 10 mM HEPES, 
pH 7.1) and incubated for 3 min at 80°C. After cooling 
down the mixture on ice, the volume was reduced by evap- 
oration at 70°C, the residue was resuspended in 4 mL of 
distilled water and centrifuged for 15 min at 13000^ and 
4°C. The supernatant was collected and the sulfate con- 
tent was then determined according to the turbidimetric 
method described by Tabatabai and Bremner [34] . 

The duplication times of the yeast cells were calcu- 
lated by fitting the equation A^o{t) = ^60o(fo) e'^' to the 
experimental data. The growth constant (kG) was esti- 
mated by expressing the growth rates (dt"^) of comple- 
mented yeasts as a function of sulfate concentrations in 
the media, and by fitting the Michaelis-Menten equation 
to the data. 

Determination of thiols, sulfate and cadmium content 

Roots and shoots were pulverized using mortar and pes- 
tle in liquid N2. Total nonprotein thiols (NPTs) and Cd 
contents were determined according to Nocito and co- 
workers [35]. Total GSH was measured according to 
Griffith [36]. 

Sulfate was extracted by homogenizing the samples in 
1:10 (w/v) ice-cold 0.1 N HNO3. After heating at SO'C for 
40 min, the extracts were filtered and the sulfate contents 
were determined according to the turbidimetric method 
described by Tabatabai and Bremner [34] . 

Sulfate influx assay and analysis of root-to-shoot 
sulfate translocation 

Sulfate influxes into the roots were measured by deter- 
mining the rates of ^^S uptake, over a 15 min pulse in 
incubation solutions labeled with the radiotracer. Briefly, 
a single plant was placed onto 400 mL of a fresh accli- 
mation solution, containing 200 \iM MgS04, supple- 
mented or not with CdCl2 at different concentrations, 
aerated and thermoregulated at 26°C. Radioactive pulses 
were started by adding ^^S -labeled Na2S04 to the uptake 
solutions. Specific activity was 4.7 kBq |imol"^. At the 
end of the pulse period, roots were excised from shoots, 
rinsed twice for 1 min in 400 mL of a 4 mM CaS04 non- 
radioactive solution at 4°C, blotted with paper towels, 
weighed, and then heated for 20 min at 80°C in 0.1 N 
HNO3 (10 mL g"^ fresh weight). Radioactivity was 
measured on aliquots of the extracting solution by li- 
quid scintillation counting in a |3 counter (LS 6000SC, 
Beckman). 
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For the analysis of root-to-shoot sulfate translocation, 
shoots were cut at 2 cm above the roots with a microtome 
blade. Xylem sap exuded from the lower cut surface was 
collected by trapping into a 1.5 mL plastic vial filled 
with a small piece of cotton for 1.5 h. The amount of 
collected sap was determined by weighing and the sul- 
fate concentration was then determined according to 
the turbidimetric method described by Tabatabai and 
Bremner [34]. 

Statistical analysis 

Statistical analysis was carried out using SigmaPlot for 
Windows version 11.0 (Systat Software, Inc.). Quantita- 
tive values are presented as mean ± standard error of the 
mean (SE). Significance values were adjusted for multiple 
comparisons using the Bonferroni correction. Statistical 
significance was at P < 0.05. Student's t-test was used to 
assess the significance of the observed differences between 
control and treated plants. Statistical significance was at 
P< 0.001. 

Results 

Cloning and functional characterization of four high-affinity 
sulfate transporter cDNAs 

Plant sulfate transporters are encoded by a multi-gene fam- 
ily whose members have specific functions in sulfate acqui- 
sition, systemic distribution and subcellular localization 
[37-39]. In this work we identified four sulfate transporter 
cDNAs expressed in B. juncea roots: one named BjSultrl;!, 
and three, with closely related sequences, named BjSul- 
trl;2a, BjSultrl;2b and BjSultrl;2c. All the cDNA-encoded 
proteins were predicted as putative high-affinity sulfate 
transporters belonging to the group 1 of the sulfate trans- 
porter family (Additional file 1). Sequence analyses re- 
vealed that the amino acid identities of these proteins with 
those of Arabidopsis belonging to the same cluster were 
86% (BjSultrl;l vs AtSultrl;l) and 94% (BjSultrl;2a/b/c vs 
AtSultrl;2), suggesting that the B. juncea and Arabidopsis 
sulfate transporters would share functions in mediating 
root sulfate uptake. In such a way BjSultrl;! could be con- 
sidered the ortholog of AtSultrl;!, whereas the three 
BjSultrl;2 cDNAs would represent three orthologous vari- 
ants of AtSultrl; 2. 

Concerning the three BjSultrl;2 forms, some additional 
data need to be taken into account. Sequence analysis 
(Additional file 2; Additional file 3) revealed that the cod- 
ing sequence of the longer variant, BjSultrl;2a, shares 
98% of nucleotide identity with Bra01S641, a gene encod- 
ing a Sultrl;2 form on the chromosome A7 of B. rapa - 
one of the two parents of B. juncea of which the genome 
has been recently sequenced [23] - and only 91% of 
nucleotide identity with Bra008340, a second form of 
Sultrl;2 found on the chromosome A2 of B. rapa. On 
the other hand the coding sequences of the shorter variants, 



BjSultrl;2b and BjSultrl;2c, share the highest identities with 
Bra008340 (95% and 99%, respectively). Unfortunately, 
we failed in finding any information about the sulfate 
transporter genes of B. nigra (the other B. juncea par- 
ent) in public genomic databases. 

The heterologous expressions of BjSultrl;!, BjSultrl;2a, 
BjSultrl;2b, and BjSultr!;2c in the yeast {Saccharomyces 
cerevisiae) double sulfate transporter mutant CP 154-7 A 
[31] were able to revert the yeast mutant phenotype, 
allowing it to grow on a minimal medium containing 
100 |iM Na2S04 as a sole sulfur source (Additional file 4), 
confirming the identity of these B. juncea clones as func- 
tional sulfate transporters. 

In order to estimate the apparent Icm for sulfate of 
each transporter we first analyzed the growth curves of 
complemented yeasts incubated in liquid media contain- 
ing different sulfate concentrations (from 0 to 100 i^M) 
as sole sulfur sources (Additional file 5). In these condi- 
tions the amount of sulfate taken up by the transporter 
and available for metabolic assimilation should be ex- 
pected to be the main limiting factor for yeast growth. 
If this were not the case - i.e. if some enzymatic activ- 
ities along the pathways of sulfate assimilation or Cys 
consumption would limit yeast growth - a gradual ac- 
cumulation of non-assimilated sulfate into the yeast 
cells should be expected. As detailed in Additional file 6, 
the sulfate content of complemented yeast cells, measured 
in the mid-log phase, did not change in the range of 
1-100 i^M sulfate external concentration, and the growth 
rate of the cells incubated in minimal media containing an 
organic sulfur source (100 |iM DL-homocysteine; HCys) 
was higher than those measured at the highest sulfate 
external concentration analyzed. Moreover, sulfate con- 
centration in the yeast cells incubated in the absence of 
sulfate was always lower than 0.05 nmol Aeoo From 
these data we can reasonably conclude that, at least in 
our conditions, the yeast growth rate is limited by sul- 
fate uptake and fits the rate of sulfate influx through the 
single heterologously expressed sulfate transporter. Thus, 
by expressing the growth rate values as a function of sul- 
fate concentrations we can calculate a growth constant, 
kc, defined as the sulfate concentration at which half of 
the maximum yeast growth rate is reached. As shown in 
Additional file 7, such a constant allows us to discriminate 
high- and low-affinity sulfate transporters, since it is 
closely related to the apparent kM for sulfate of the trans- 
porters. Least square fittings (Figure 1) revealed that the 
growing isotherms of the four complemented yeasts can 
be properly described by single hyperbolic Michaelis- 
Menten functions, with kc for sulfate in the micromolar 
range, indicating that these proteins are high-affinity sul- 
fate transporters; in particular the kg values were 5.46 ± 
0.22 \iM (BjSultrl;l), 1.74 ±0.05 |iM (BjSultrl;2a), 1.73 ± 
0.07 |iM (BjSultrl;2b), and 1.74 ± 0.05 ^iM (BjSultrl;2c). 
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Figure 1 Estimation of the growth constant (kc) dependent on sulfate of the yeasts expressing the Brassica juncea sulfate transporters. 

The duplication times (dt) of tlie complemented yeast cells were calculated by fitting the equation A(a>(t) = /\6oo(fo) e^' to the experimental data 
reported in Additional file 5. kg was estimated by expressing the growth rates (dt^^) of complemented yeasts as a function of sulfate concentrations 
in the media, and by fitting the Michaelis-Menten equation to the data. Data points and error bars are means and SE of two experiments run in 
triplicate (n = 6). 



Effect of Cd exposure and sulfate limitation on sulfate 
uptake and sulfur allocation in Brassica juncea plants 

All the data presented in this paragraph derived from ex- 
periments aimed at comparing environmental conditions 
(Cd exposure and sulfate limitation) in which sulfate up- 
take induction should occur. For these purposes plants 
were exposed to 10 and 25 \iM Cd^^ for 2 days or grown 
under sulfate limitation (50 and 10 \iM SO4 ^ ) for a 10- 
day period. Control plants were grown at 200 SO4 
in the absence of Cd. 

Cadmium exposure neither significantly influenced 
the growth of shoots and roots, nor produced any ap- 
parent symptom of stress; conversely, lowering sulfate 
concentration in the growing solution significantly in- 
creased root growth without affecting shoots, as indi- 
cated by the values of the shoot/ root ratio which 
decreased from 3.82 to 2.20 (Table 1). The total amount 
of Cd retained by roots increased as the metal concen- 
tration in the external medium did, whereas it reached 
similar values in the shoots of plants grown at 10 and 
25 mM Cd^+ (Table 1). 



Cadmium exposure and sulfate limitation deeply af- 
fected the sulfate uptake capacity of the root, as indi- 
cated by the values of ^^S -sulfate uptake, measured at 
200 \iM S04^" external concentration (Figure 2A, B). In 
Cd exposed plants the rate of sulfate uptake increased 
up to 0.9-fold with respect to the untreated control at 
the highest Cd external concentration (25 mM). Similar 
behaviors were observed in sulfur- starved plants, in which 
the rate of sulfate uptake increased as the sulfate con- 
centration in the external medium decreased, reaching 
values 1.2-fold higher than in sulfur-sufficient control 
(200 \iM S04^ ). These trends were closely associated 
to changes in the transcript level of BjSultrl;! and in 
the cumulative amount of the three transcripts of the 
BjSultrl;2 variants {BjSultrl;2 pool), which significantly 
accumulated as the severity of the stresses increased 
(Figure 2C, D). 

Taken as a whole these preliminary results indicate that 
48 h Cd exposure and 10-day sulfate limitation produced 
similar induction of sulfate uptake. Since such effects 
should presumably be related to changes in the sulfur 
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Table 1 Dry weight of roots and shoots and Cd accumulation 


Experimental condition 


Dry weight 


Shoots/Roots 


Cd^* content 




Roots 


Shoots 




Roots 


Shoots 




g 






\xmo\ g~ 


DW 


Control 


0.063 ± 0.003 (a) 


0.241 ±0.010 (a) 


3.82 


ND 


ND 


10 mM Cd^* 


0.069 ± 0.003 (a) 


0.252 ±0.011 (a) 


3.65 


25.81 ±1.18 (a) 


5.66 ±0.25 (a) 


25 mM Cd^+ 


0.066 ± 0.004 (a) 


0.231 ±0.012 (a) 


3.50 


97.33 ± 3.99 (b) 


5.00 ±0.22 (a) 


50 |jM S04^" 


0.1 12 ±0.005 (b) 


0.251 ±0.011 (a) 


2.24 


ND 


ND 


10 |jM S04^" 


0.110±0.005 (b) 


0.243 ±0.011 (a) 


2.20 


ND 


ND 



Plants were exposed to different Cd concentrations (10 and 25 pM) for 48 h or grown under different sulfate concentrations (50 and 10 pM) for 10 days. Control 
plants were grown under 200 S04^" and were not exposed to Cd. Cadmium content was measured by ICP-IVIS. Values are means ±SE of three experiments 
run in triplicate (n = 9). Different letters indicate significant differences (P<0.05). ND, not detectable. 



nutritional status of the plants, we analyzed the levels of 
NPTs, GSH and sulfate of both roots and shoots, as- 
suming the pools of these intermediates as the main 
diagnostic indicators of the sulfur nutritional status. 

Cd exposure produced significant changes in the NPT 
levels of the root, which progressively increased as Cd 
concentration in the external medium did (Figure 3A), 
whilst at the same time a decrease of the total GSH 
pools was observed (about 30% with respect to the con- 
trol in all analyzed conditions; Figure 3B). Such a trend 
was probably related to PC biosynthesis and accumula- 
tion according to the progressive increase in Cd root 



content (Table 1). The sulfate pools of the root were not 
affected by Cd exposure (Figure 3C). Quite similar be- 
haviors were observed in the shoots of Cd exposed 
plants, since the NPT levels increased with Cd concen- 
tration in the external medium and the sulfate concen- 
tration was not affected by the presence of the metal; 
however, a Cd-dependent increase in the GSH levels was 
observed (Figure 3 A, B, C). 

As expected, a stepwise contraction in the levels of all 
the diagnostic indicators was observed in the root of 
plants grown for 10 days under sulfate limitation. In- 
deed, NPT, GSH and sulfate contents measured in the 
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Figure 2 Changes in sulfate uptake capacity of Brassica juncea roots. Plants were exposed to different Cd concentrations for 48 h (A, C) or 

grown under different sulfate concentrations for 10 days (B, D). (A, B) Sulfate influxes were evaluated by measuring the rate of ^^S04^" absorption 
Into roots of Intact plants over a 15 min pulse. The Incubation solutions contained 200 ^iM S04^". Bars and error bars are means and SE of three 
experiments run In triplicate (n = 9). Different letters Indicate significant differences (P<0.05). (C, D) Seml-quantltatlve RT-PCR analysis of 
BJSultrl;! and BjSultrl:2 gene expression. PCRs were carried out for 24 cycles where cDNAs were exponentially amplified. For BjSuitrl;2 pool, 
primers were designed on conserved sequences of the three BJSultrl:2 variants, and gave overlapping amplification products of 1046 bp. PCR 
products were separated in agarose gel and stained with SYBR Green I. Signals were detected using a laser scanner with 532 nm laser and 
526 nm filter. BjTub, tubulin. A representative set of data from three independent experiments Is given. 



Lancilli ef al. BMC Plant Biology 2014, 14:132 
http://www.biomedcentral.com/1471-2229/14/132 



Page 8 of 15 



5000 
4000 
3000 
§ 2000 
1000 

IT 

w 0 

o 

o 
E 
c 

I 1000 
c 

800 



b 

I 



10 

[ Cd'" ] (MM) 



25 



600 



400 



200 



a 

LiLl 



200 



50 

, 2- 



10 



[ so/- ] (MM) 




I 

w 



c 




25 


1 




20 






15 






10 




? 






LL 


5 




O) 






O 


0 



800 



600 



400 



200 



0 10 25 
[ Cd'" ] (MM) 



a 



200 



50 



10 



[SO/-] (MM) 




0 10 25 



[ Cd" ] (UM) 




[ SO/- ] (MM) 



Figure 3 Effects of Cd exposure and sulfate limitation on the sulfur nutritional status of Brassica juncea plants. Plants were exposed to 
different Cd concentrations for 48 h (A, B, C) or grown under different sulfate concentrations for 10 days (D, E, F). (A, D) NPT contents of roots 
(black bars) and shoots (grey bars) are expressed as GSH equivalents. (B, E) Total GSH contents of roots (black bars) and shoots (grey bars). 
(C, F) Sulfate contents of roots (black bars) and shoots (grey bars). Bars and error bars are means and SE of three experiments run in triplicate 
(n = 9). Different letters indicate significant differences between treatments (P< 0.05). ND, not detectable. 



root tissues dramatically decreased as sulfate availability 
in the external medium did (Figure 3D, E, F). Following 
sulfate limitation, sulfate content of the shoot steadily 
decreased, reaching the minimal value at 10 mM S04^" 
external concentration; differently, NPT and GSH levels 
did not significantly change when we lowered sulfate ex- 
ternal concentration from 200 to 50 ^M, whilst a sharp 
decrease in the level of these compounds was observed 
by moving toward the lowest (10 mM) sulfate concentra- 
tion analyzed (Figure 3D, E, F). 

We also analyzed the dynamic of root-to-shoot sulfate 
translocation by measuring the concentration of the anion 
in the xylem sap of Cd-exposed or sulfur-starved plants. 
In these experiments, sulfate translocation was estimated 
as the amount of sulfate ions loaded and transported in 
the xylem sap for 1.5 h. Results indicate that the amount 
of sulfate ions transported in the xylem sap progressively 
increased following Cd exposure (Figure 4A); differently, 
sulfate translocation increased when shifting sulfate exter- 
nal concentration from 200 to 50 |iM, and sharply de- 
creased when moving toward the lowest (10 [iM) sulfate 
concentration analyzed (Figure 4B). 



Quantitative analysis of the expression of the three 
BjSultrl;2 variants 

Since the three BjSultrl;2 forms are not polymorphic 
enough to be distinguished by means of a simple PGR 
(Additional file 8), we developed a suitable method to study 
changes in their expression by coupling semi-quantitative 
RT-PGR analysis with the use of an opportune restriction 
enzyme. 

Sequence analysis revealed that the three BJSultrl;2 
cDNAs have restriction site polymorphisms for the Glal 
endonuclease, which enabled us to discriminate the 
three variants after digestion. As detailed in Figure 5A: 
i) BjSultrl;2a (1968 bp) is not cut by Glal; ii) BjSultrl;2b 
(1959 bp) is cut by Glal 1752 bp downstream of the start 
codon; iii) BjSultrl;2c (1959 bp) is cut twice by Glal, 
1098 and 1752 bp downstream of the start codon. As a 
consequence the digestion of the cDNA clones with Glal 
produces characteristic restriction patterns with some 
diagnostic bands useful to discriminate the three forms 
(Figure 5B). The characteristic undigested 1968 bp band 
is a diagnostic marker of BjSultrl;2a presence, the 1752 bp 
fragment only results from the digestion of BjSultrl;2b, 
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Figure 4 Effects of Cd exposure and sulfate limitation on sulfate 
translocation in Brassica juncea plants. Plants were exposed to 
different Cd concentrations for 48 h (A) or grown under different sulfate 
concentrations for 10 days (B). At the end of the experimental periods, 
shoots were separated from roots and the xylem sap exuded from the 
cut (root side) surface was collected to be analyzed for sulfate content. 
Bars and error bars are means and SE of three experiments run in 
triplicate (n = 9). Different letters indicate significant differences (P < 0.05). 



whilst both the 1098 and 654 bp bands are specifically 
produced following the digestion of BjSultrl;2c. Finally 
the 207 bp band is a digestion product shared among 
BjSultrl;2b and BjSultrl;2c, and therefore does not give 
any result useful for our purposes. 

Starting from this rationale, we designed a couple of 
primers amplifying at the same time the entire open read- 
ing frames of the three clones with the same efficiency (data 
not shown), and we used these oligos for the semi- 
quantitative RT-PCR analysis of the effects of 48 h exposure 



to 25 nM Cd^"^ or 10-day sulfate limitation (10 |iM S04^ ) 
on the expression of the three BjSultrl;2 variants. A restric- 
tion analysis using the Clal endonuclease followed the amp- 
lification reactions. 

Results show that the cumulative amount of the 
BjSultrl;2 transcripts in the roots was significantly 
higher in Cd exposed plants than in the untreated control 
ones (+217%, Figure 6A) as already shown in Figure 2C. 
Such a behavior resulted from changes in the expression 
of BjSultrl;2b and BjSultrl;2c only (Figure 6A). In fact, 
the densitometric analysis of each diagnostic band indi- 
cated that BJSultrl;2b and BjSultrl;2c transcript levels sig- 
nificantly increased by 585% and 301%, respectively, whilst 
the BjSultrl;2a expression was not affected by Cd expos- 
ure (Figure 6B, Additional file 9). Similar behaviors were 
observed by analyzing changes in the expression pattern 
of the three BjSultrl;2 forms in plants exposed for 48 h to 
a lower (10 \iM) Cd concentration (Additional file 10). In 
fact, also in this condition the response to Cd, though to a 
lesser extent, was only ascribable to specific increases in 
the relative amount of BjSultrl;2b (+150%) and BjSultrl;2c 
(+75%) transcript. 

By contrast, in the case of sulfate limitation, the increase 
in the cumulative amount of the BjSultrl;2 transcripts 
(+455% with respect to the sulfur-sufficient control) was 
ascribable to changes in the transcript levels of all the 
three forms (Figure 6C). In particular, the BjSultrl,2a, 
BjSultrl;2b, and BJSultrl;2c transcript levels signifi- 
cantly increased by 371%, 483%, and 618%, respectively 
(Figure 6D, Additional file 9). 

Discussion 

Brassica juncea (L.) Czern & Coss (AABB, n = 18) is 
believed to have originated from the interspecific hy- 
bridization of two base "diploid" genomes provided by 
Brassica rapa L. (AA; n = 10) and Brassica nigra L. (BB; 
« = 8) [24,40]. Both the diploid parents are thought to 
be ancient polyploids since they still exhibit highly repli- 
cated genomes, each containing three paralogous subge- 
nomes closely related to that of Arabidopsis thaliana 
[21-23]. In spite of the whole-genome triplication event - 
thought to have occurred between 13 and 17 million years 
ago - most comparative studies have shown that the num- 
ber of each ancestral gene retained in the genome of the 
modern diploid Brassica is variable, since paralogous 
regions exhibit interspersed gene losses and insertions. 
Interestingly, in the recently sequenced B. rapa genome 
the extent of gene loss among triplicated genome seg- 
ments varies, with one of the three copies consistently 
retaining a disproportionately large fraction of the genes 
expected to have been present in its ancestor [23]. Such 
evolutionary events are thought to be the biological basis 
of the immense plasticity of Brassica species and may have 
led to a diversification of the genes retained in more than 
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Figure 5 Restriction analysis of the three BjSultr1;2 cDNAs. (A) The three BjSultr1;2 variants have restriction site polymorphisms for the Cla I 
endonuclease. Black arrows indicate the relative position of Cla I restriction sites in each open reading frame. The expected lengths of the restriction 
fragments obtained after digestion with Clal are indicated. (B) Charaaeristic restriction patterns obtained from the digestion of each cDNA with Clal. 
Single cDNAs were obtained by PCR using a recombinant plasmid, containing a unique Bj5ultrl:2 clone, as template, u, undigested; d, digested. 



one copy, in terms of function and/or expression. Search- 
ing for orthologs of the Arabidopsis high-affinity sulfate 
transporter genes involved in sulfate uptake and retained 
in the genome of B. rapa - one of the two parents of 
B. juncea - revealed the existence of three distinct loci, 
annotated as Bra022623, Bm01S641 and Bm008340. The 
first locus encodes for a putative high-affinity sulfate trans- 
porter closely related to Arabidopsis AtSultrl;!, whilst the 
other two loci encode for two different forms of a high- 
affinity sulfate transporter functional related to Arabidopsis 
AtSultrl;2, indicating these gene loci as paralogs. As 
expected, a much more complex picture was found in 
the allopolyploid B. juncea in which we were able to 
identify an orthologous form of AtSultrl;! (BJSultrl;!) 
and three orthologous forms of AtSultrl;2 {BjSultrl;2a/b/c). 
From the results obtained by the sequence analysis, and 
in the absence of any other information so far available 
about the B. nigra genome, we can reasonably suppose: 
i) BjSultrl;2a as the ortholog of Bra01S641 on the gen- 
ome A or B of B. juncea; ii) BjSultrl;2blc as allelic variants 
orthologous of Bra008340 on the A or B genome of B. 
juncea or a homeologous gene pair related to Bra008340 
on the A and B genomes of B. juncea. Moreover, since 
the progeny derived from self-fertilization inherited all 
the three BjSultrl;2 variants (data not shown), it seems 



likely to exclude that BjSultrl;2b and BjSultrl;2c would 
be allelic, making plausible the hypothesis they are instead 
present at different homeologous gene loci on A and B ge- 
nomes, each in homozygous configuration; otherwise, a 
simple mendelian segregation would be observed. In any 
case, since the three BjSultrl;2 forms would share a com- 
mon ancestor gene, they may either have retained their 
original functions and expressions, or - as it is often the 
case - have accumulated deleterious mutations or have 
evolved novel gene interactions through the processes of 
sub-functionalization and/or neo-functionalization [25,26]. 

Results of complementation tests in the yeast mutant 
strain CP154-7A, defective in its two sulfate transporters 
and thus unable to grow on media containing low concen- 
trations of sulfate as the sole sulfur source [31], proved 
the capacity of BjSultrl;! and BjSultrl;2a/b/c to transport 
sulfate ions across the plasma-membrane (Additional 
file 4). Moreover, kinetic analysis of the growth (G) iso- 
therms of complemented yeasts, revealed that BjSultrl;! 
and BjSultrl;2a/b/c have high affinities for sulfate, as re- 
vealed by the 1<.q values similar to the apparent Icm values 
of other plant high-affinity sulfate transporters [9,41-44] 
indicating that all the B. juncea clones have retained 
their functions. It is also worthy of note that the sulfate 
transporter BjSultrl;l has an apparent affinity for sulfate 
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Figure 6 RT-PCR analyses of the three BjSultr1;2 forms in the roots of Brassica juncea grown under Cd exposure or sulfate limitation. 

Plants were exposed to 25 ^iIVl Cd^^" for 48 h (A, B) or grown under 10 ^iM S04^" for 10 days (C, D). (A, C) The entire ORFs of the three BJSultrl:2 
forms were amplified and PCR products, digested (d) or not digested (u) with Clal endonuclease, were electrophoresed on agarose gel and 
stained with SYBR Green I. cDNA loading was normalized using BJTub as an internal control. Signals were detected using a laser scanner with 
532 nm laser and 526 nm filter. (B, D) Densitometric analysis. Arrows indicate the relative position of each electrophoretic band obtained after 
digestion of PCR products with Clal. A representative set of data from three independent experiments is given. For statistical analysis, see 
Additional file 9. 



(kc3 = 5.46 |iM) three times lower than those of the three 
BjSultrl;2 forms (Figure 1). Finally, phylogenetic ana- 
lyses indicate these transporters as functionally related 
to AtSultrl;! and AtSultrl;2 and then we can infer their 
probable function in mediating root sulfate uptake from 
the soil solution [14,43,44]. 

Physiological analysis reveals that Cd exposure as well 
as sulfate limitation induces sulfate uptake in B. juncea 
roots. Such behaviors seem to be related to the induc- 
tion of high-affinity sulfate transporters belonging to the 
group I, as indicated by the increase in the transcript 
levels of BjSultrl;! and BjSultrl;2 pool (Figure 2). Thus - 
at the physiological level - B. juncea also retains the typ- 
ical responses of sulfate uptalce to Cd or sulfur shortage 
[9,15]. The apparent quantitative discrepancy between the 
changes in the transcript levels of BjSultrl;! and BjSultrl;2 
pool and the resulting increases in sulfate uptake may be 
due to additional regulatory mechanisms working in paral- 
lel with the transcriptional control of the high- affinity sul- 
fate transporter genes [45]. 

Considering the current model of demand-driven regu- 
lation of sulfate uptake, such inductions should be related 
to the sulfur nutritional status reached by plants in 
the two growing conditions. Changes in the amounts 
of sulfur-containing compounds that we assume as the 



main diagnostic indicators of the sulfur nutritional status 
of root and shoot clearly indicate that Cd exposure and 
sulfate limitation influence sulfur allocation throughout 
the whole plant, generating deeply different local nutri- 
tional states which make it difficult to individuate an un- 
equivocal and common nutritional signal related to the 
expression of sulfate transporters (Figure 3). Indeed, as 
expected, lowering sulfate concentration in the external 
medium necessarily results in a significant contraction 
of all the analyzed sulfur pools of the roots, whilst Cd 
stress produces a typical increase in the level of NPTs, 
probably due to the activation of GSH-dependent PC 
biosynthesis, without affecting the sulfate content of 
the roots. 

Root responses to Cd exposure seem to be due to 
homeostatic mechanisms driven by increases in the sul- 
fur need of the plants, since, as previously reported [9], 
the effect of Cd on sulfate uptake capacity is closely re- 
lated to the NTP levels of both root and shoot and then 
to the strength of the Cd-induced additional sink for thiols 
(Additional file 11). Under Cd exposure the NPT levels in 
plant tissues significantly increase, reaching values 8.5 
(root) and 1.3 (shoot) fold higher than in the control at 
the highest concentration analyzed. Since the sulfate pools 
of both root and shoot seem not to be affected by Cd 
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exposure, it appears clear that the additional sulfur re- 
quired to sustain thiol biosynthesis necessarily derives from 
the activation of sulfate uptake. On the other hand, root re- 
sponses to sulfate shortage appear likely to be dependent 
on the need for allocating the limiting nutrient in the best 
and most efficient way. In fact, the stepwise decrease in the 
sulfate content of the roots seems to be related not only to 
a decrease in sulfate availability, but also to a transient acti- 
vation of sulfate translocation making shoots less sensitive 
to sulfate deficiency (Figure 4B) [46] . Noticeably, Cd expos- 
ure neither decreases root sulfate content (Figure 3C) nor 
inhibits root sulfate uptake (Figure 2A) and sulfate trans- 
porter gene expression (Figure 2C), but rather significantly 
enhances sulfate translocation (Figure 4A) and its metabol- 
ism as shown by the significant increases in the GSH levels 
of the shoot (Figure 3B). Such an effect should be related 
to the activation of both PC biosynthesis and mechanisms 
involved in controlling oxidative damage due to the accu- 
mulation of free Cd ions in the leaves [47,48]. Moreover, 
from these results we can also speculate that the over- 
accumulation of GSH in the shoot could help roots in de- 
toxifying Cd through reallocating mechanisms involving 
phloem translocation, as previously reported in other spe- 
cies [49,50]. In this way, the excess of sulfate taken up by 
roots would partiy bypass root assimilation to be directly 
metabolized in organs less affected by Cd stress, without 
however affecting the root sulfate pool (Figure 3C). 

Taken as a whole our data clearly show that dissimilar 
nutritional and metabolic states may result in quite similar 
responses in sulfate uptake, suggesting that multi-signalling 
pathways may control the expression of the high-affinity 
sulfate transporters of the roots. Moreover, the fact that the 
negative relationships between the levels of nutritional sig- 
nals (sulfate and GSH) and sulfate uptake capacity of the 
roots, existing in sulfur-starved plants, are not found in the 
Cd exposed ones, seems to further support this conclusion. 

Although it is difficult to indicate unambiguous nutri- 
tional signals, we can make some educated guesses on 
the general structure of the hypothetical signaling path- 
ways involved in the modulation of sulfate uptake by 
analyzing the expression pattern of the three BjSultrl;2 
forms under Cd exposure and sulfate limitation. Since 
BJSultrl;2a seems to have lost its capacity to respond to 
Cd stress, but, at the same time, retains its response to 
sulfate shortage, we can speculate about the existence of 
at least two distinct signal transduction pathways. The 
first one modulates root sulfate uptake, ensuring ad- 
equate nutrient supply when plants experience lowering 
in the sulfate concentration of the soil solution, probably 
through a c«-acting sulfur responsive element as previ- 
ously suggested [19], whilst the second one we postulate 
is likely to be involved in meeting sulfate uptake with 
the plant metabolic sulfur demand, which may increase 
following heavy-metal stress. All the three BjSultrl;2 forms 



are controlled by the first pathway as indicated by the ana- 
lysis of their respective contribution to the increase in the 
cumulative amount of the BjSultrl;2 transcripts under sul- 
fur starvation, but only two forms {BJSultrl;2b/c) seem to 
have retained the ancient characteristic to be controlled by 
the second regulatory pathway. In this context the differen- 
tial transcriptional behaviors of the three BjSultrl;2 forms 
could be explained by hypothesizing the presence of both 
"cadmium-" and "sulfur-sensitive" regions in the promoter 
of an ancient Brassica Sultrl;2 form, which - following 
polyploidization - may have evolved in sub-functionalized 
forms, whose combined actions result in molecular and 
physiological responses to Cd exposure and sulfate limita- 
tion similar to those known in species with non-redundant 
genome [8,9,15]. If this were not the case an interference 
of the metal with the signal transduction pathways in- 
volved in the regulation of sulfate uptake should be pos- 
tulated, as suggested in the recent paper of Shahbaz and 
co-workers [51], in which they extensively discuss the 
effect of copper accumulation on sulfur metabolism-related 
gene expression. However, copper stress in Brassica seems 
to produce significant increases in both sulfate uptake and 
tissue sulfate content without substantially altering plant 
sulfur demand [52], differentiy from Cd exposure which 
produces increases in sulfate uptake closely related to 
the strength of the additional sink for thiols it induces 
(Additional file 11). Moreover, if any sort of Cd interfer- 
ence occurs we should also suppose the existence of at 
least two signal transduction pathways controlling the 
expression of BjSultrl;2a/b/c under Cd exposure: the 
first inhibited by Cd, and the second Cd insensitive. Fi- 
nally, we cannot exclude that other molecular mechanisms 
may be involved in the differential expression of the three 
BjSultrl;2 forms under Cd stress, as for instance those 
suggested for rice phosphate transporters [53]; it could be 
interesting to investigate if the short 9 bp-insertion lo- 
cated at the 5' end of the BjSultrl;2a coding sequence 
(Additional file 8) can play a role in the regulation of its 
expression. 

Conclusions 

Taken as a whole our data agree with the main molecu- 
lar and physiological evidence obtained in Arabidopsis 
which support the idea that the regulation of AtSultrljl 
and AtSultrl;2 - the two transporters mediating sulfate 
uptake from the soil solution - must necessarily involve 
independent signaling pathways, as extensively shown by 
Rouached and co-workers [15,54]. Moreover, we can also 
conclude that different sulfur nutritional and metabolic 
conditions may be perceived by a single sulfate trans- 
porter gene. Such a finding reveals that the mechanisms 
involved in sulfate uptake regulation may be more com- 
plex than previously thought, and partially accounts for 
the lack of unambiguous nutritional signals, since the 
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activity of each transporter may result from a complex 
interplay among multiple regulatory pathways. 

Additional files 



Additional file 1: Dendrogram showing sulfate transporter family 
of Arabidopsis thaliana and high-affinity sulfate transporters of 
Brassica juncea. The dendrogram was constructed on the bases of 
amino acid sequences using MEGA 5.05 software. Accession numbers for 
A. thaliana (TAIR; http://www.arabidopsis.org/) are: AtSultrl:!, At4g08620; 
AtSultrl:2, Atlg78000: AtSultrl:3, Atlg22150: AtSultr2:l, At5gl0180: 
AtSultr2:2, Atlg77990: AtSultr3:l, At3g51895: AtSultr3:2 At4g02700: 
AtSultr3:3, Atlg23090: AtSultr3;4, At3gl5990: AtSultr3:5, At5gl9600: 
AtSultr4:l, At5gl3550: AtSultr4:2, At3gl2520. Accession numbers for B. 
juncea (GenBank; http://www.ncbi.nlm.nih.gov/genbank/) are: BjSultrl;!, 
JX896426; BjSultrl ;2a, JX896427; BjSultrl ;2b, JX896428: BjSultrl ;2c, 
JX896429. 

Additional file 2: Nucleotide identity (%) between the coding 
sequences of the three BjSuitrl;2 variants and other Brassica 
Sultrl;2 coding sequences. 

Additional file 3: Dendrogram showing high affinity sulfate 
transporters of Arabidopsis thaliana, Brassica juncea, Brassica napus, 
and Brassica rapa. The dendrogram was constructed on the bases of 
amino acid sequences using MEGA 5.05 software. Accession numbers for 
A. thaliana (TAIR; httpy/www.arabidopsis.org/) are: AtSultrl;!, At4g08620; 
AtSultrl ;2, Atlg78000. Accession numbers for B. juncea and B. napus 
(GenBank; http://www.ncbi.nlm.nih.gov/genbank/) are: BjSultrl;], JX896426; 
BjSultrl ;2a, JX896427; BjSultrl ;2b, JX896428; BjSultrl ;2c, JX896429; BnSultrl ;1 , 
AJ415460; BnSultrl ;2, AJ31 1388. Accession numbers for ft rapa (BRAD; 
http://brassicadb.org/brad/) are: Bra022623; Bra015641; Bra008340. 

Additional file 4: Phenotypic complementation of the yeast double 
sulfate transporter mutant CP154-7A by the sulfate transporters of 
Brassica juncea. Yeast mutant cells expressing BjSultrl;], BjSultrl:2a, 
BjSultrl:2b, and BjSultrl:2c under the control of the galactose-inducible 
GALIO promoter or harboring the empty pESC-TRP vector were grown at 
28°C for 3 d on a minus-sulfur minimal medium (-S) or on minimal 
media containing 100 (jM sulfate (SO^^^) or 100 (jM DL-homocysteine 
(HCys) as sole sulfur sources. 

Additional file 5: Growth curves of complemented yeast cells. (A, B, 

C, D) Complemented yeasts were incubated at 28°C for 25 h in liquid 
media containing different sulfate concentrations (• 0 (jM; o 1 (jM; T 
2.5 mM; a 5 tjM; ■ 7.5 |jM; □ 10 ^M; ^ 25 mM; O 50 pM; ▲ 100 pM) or 

100 (jM HCys (V) as sole sulfur source. Absorbance was measured at 
600 nm (A^o) along time. Data points and error bars are means and SE 
of two experiments performed in triplicate (n = 6). 

Additional file 6: Sulfate content in complemented yeast cells. 

Complemented yeast cells expressing BjSultrl;l, BjSultrl ;2a, BjSultrl ;2b, 
and BjSultrl ;2c were incubated in liquid media containing different 
sulfate concentrations as sole sulfur source. At the end of the incubation 
period yeasts were harvested and processed to determine their sulfate 
content. Values are means + SE of two experiments run in triplicate (n = 6). 

Additional file 7: Growth analysis of CP154-7A yeast mutant 
expressing an high- or a low-affinity sulfate transporter of Zea mays 
(ZmSTl;!) or Arabidopsis thaliana (AtSultr2;l), respectively. ZmSTl;! 
and AtSultr2;l coding sequences were amplified by RT-PCR from total 
RNA isolated from maize and Arabidopsis roots, respectively, and cloned 
in the pESC-TRP vector as described in Methods. (A) Complemented yeast 
cells were incubated in liquid media containing two sulfate concentrations 
(0.1 or 0.5 mM) as sole sulfur source. Yeast growth was monitored by 
measuring the/\6oo nm at different times. (•) ZmSTl ;1 at 0.1 mM sulfate; 
(o) ZmSTl ;1 at 05 mM sulfate; (T) AtSultr2;l at 0.1 mM sulfate; (V) 
AtSultr2;l at 0.5 mM sulfate. (B) Growth curves of yeast cells expressing 
ZmSTl ;1 (• 0 IJM; o 1 |jM; T 2.5 ^M; A 5 mM; ■ 7.5 (jM; □ 1 0 (jM; ♦ 25 (jM; 
O 50 [iM; ▲ 100 |jM). (C) Growth curves of yeast cells expressing AtSultr2;l 
(• 0 mM; o 0.05 mM; T 0.1 mM; A 0.1 5 mM; ■ 02 mM; □ 025 mM; ♦ 05 
mM; O 1 mM; ▲ 1 5 mM; V 2 mM; # 2.5 mM; O 3 mM). (D, E) Estimation 
of the grovrth constant (kcjj for sulfate. The duplication times (dt) of the 



complemented yeast cells were calculated by fitting the equation /Igooft) = 
'^6oo(to) to the experimental data reported in B and C. ks was determined 
by expressing the growth rates (dt"') of complemented yeasts as a function 
of sulfate concentrations in the media, and by fitting the Michaelis-Menten 
equation to the data. Results reveal that the k^ values for sulfate were similar 
to the kivi values measured for each sulfate transporters by using conventional 
methods (Nocito etal Plant Physiol, 2006 141:1 138-1 14S; Takahashi ef al. Plant 
J, 2000 23:1 71-182). Data points and error bars are means and SE of two 
experiments performed in triplicate (n = 6). 

Additional file 8: Alignment of nucleotide sequences of the three 
BjSultrl;2 forms. Shared nucleotides are highlighted in grey. 

Additional file 9: Changes in the transcript relative amount of the 
three BjSultrl;2 forms in the roots of Brassica juncea grown under 
Cd exposure or sulfate limitation. Plants were exposed to 25 pM Cd^* 
for 48 h (+Cd) or grown under 10 pM S04^" for 10 days (-S). Control 
plants were grown under 200 pM S04^" and were not exposed to 
cadmium. The entire OREs of the three BjSultrl:2 forms were amplified 
and PCR products were digested with Clal endonuclease, electrophoresed on 
agarose gel, and finally stained with SYBR Green I. Signals were detected using 
a laser scanner with 532 nm laser and 526 nm filter and densitometrically 
analyzed using ImageJ 1.46 software. cDNA loading was normalized using 
BjTub as an internal control. Bars and error bars are means and SE of three 
independent experiments run in triplicate (n = 9). Asterisks indicate significant 
differences between control and treated plants {P< 0.001). 

Additional file 10: RT-PCR analyses of the three BjSultrl ;2 forms in 
the roots of Brassica juncea exposed to 10 pM Cd. Plants were 
exposed or not to 10 pM Cd^* for 48 h. (A) The entire OREs of the three 
BJ5u!trl;2 forms were amplified and PCR products were digested with Clal 
endonuclease, electrophoresed on agarose gel, and finally stained with 
SYBR Green I. cDNA loading was normalized using BjTub as an internal 
control. Signals were detected using a laser scanner with 532 nm laser 
and 526 nm filter. A representative set of data from three independent 
experiments is given. (B) Densitometric analysis. Arrows indicate the 
relative position of each electrophoretic band obtained after digestion of 
PCR products with Clal. (C) Statistical analysis. Bars and error bars are 
means and SE of three independent experiments run in triplicate (n = 9). 
Asterisks indicate significant differences between control and treated 
plants {P<a001). 

Additional file 11: Relationship between NPT content and sulfate 
uptake capacity in plant of Brassica juncea exposed to different Cd 
concentrations. Plants were exposed for 48 h to differentCd^* 
concentrations: 0 (white), 10 (grey), and 25 (black) pM. Circles, roots; 
triangles, shoots. Data points and error bars are means and SE of three 
experiments run in triplicate (n = 9). 
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